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Diet is central for understanding hominin evolution, adaptation, and 
environmental exploitation, but Paleolithic plant remains are scarce. 
A unique macrobotanical assemblage of 55 food plant taxa from the 
Acheulian site of Gesher Benot Ya'aqov, Israel includes seeds, fruits, 
nuts, vegetables, and plants producing underground storage or- 
gans. The food plant remains were part of a diet that also included 
aquatic and terrestrial fauna. This diverse assemblage, 780,000 y 
old, reflects a varied plant diet, staple plant foods, environmental 
knowledge, seasonality, and the use of fire in food processing. It 
provides insight into the wide spectrum of the diet of mid-Pleisto- 
cene hominins, enhancing our understanding of their adaptation 
from the perspective of subsistence. Our results shed light on hom- 
inin abilities to adjust to new environments, facilitating population 
diffusion and colonization beyond Africa. We reconstruct the major 
vegetal foodstuffs, while considering the possibility of some detox- 
ification by fire. The site, located in the Levantine Corridor through 
which several hominin waves dispersed out of Africa, provides a 
unique opportunity to study mid-Pleistocene vegetal diet and is 
crucial for understanding subsistence aspects of hominin dispersal 
and the transition from an African-based to a Eurasian diet. 
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D iet is central for understanding the evolution, adaptation, en- 
vironmental exploitation, cognition, technology, and survival of 
prehistoric hominins. Reconstructions of Acheulian diets are based 
on skeletal material (1), isotopic signatures (2), ecological models 
reconstructing African paleoenvironments (3), comparative studies 
of primate behavior, especially that of chimpanzees and bonobos 
(4), and the diets of modern hunter-gatherers (5). 

Direct data on Paleolithic plant diets are scarce, because plant 
remains are perishable, and most information is circumstantial 
(e.g., isotopic ratios reflecting C3/C4 plant taxa use relations) and 
insufficient for detailed reconstruction (1). Direct but limited evi- 
dence for plant consumption is sometimes found, however, in 
calculus (6). Earlier overemphasis of animal proteins and fats in 
reconstructions of prehistoric diet because of their better-preserved 
remains has been somewhat moderated by using ethnographic 
analogies (7). Recently, marine resources have also been consid- 
ered (8) but are relevant only for coastal or aquatic-related sites. 
Overall, understanding the role of plants in early hominin diets has 
been based on meager direct evidence from only a few sites, such 
as Kalambo Falls in Africa (9) and Karlich in Europe (10). 

The waterlogged Acheulian site of Gesher Benot Ya'aqov (GBY) 
yielded many well-preserved macrobotanical remains, including 
wood and bark (11); fruits, nuts, and seeds (12, 13); and pollen (14). 
Many inedible food plant seeds belong to species that have other 
plant food organs, such as vegetables and underground storage or- 
gans (USOs) (13, 15). These plant assemblages, originating in the 
Levantine Corridor through which several hominin waves dispersed 
out of Africa, create a special opportunity to study hominin vegetal 
diet during Early-Middle Pleistocene times. GBY is, thus, a key site 
for understanding hominin dispersal and colonization out of Africa 
from the perspective of plant food gathering, nutrition, and sea- 
sonality, illuminating the transition from an African-based diet to a 


Eurasian one. Our focus on plant foods stems from GBY’s excep- 
tionally rich botanical remains as well as comparisons with diets of 
current hunter-gatherers (16) and wild food plant gathering in tra- 
ditional Near Eastern societies. 

The Hula Valley (196 km 2 ) (17) and the surrounding moun- 
tains form a catchment area of ~1,500 km 2 , harboring over 300 
food plant species (18) (Table SI). The immediate environ- 
mental setting of GBY included three different habitats: ( i ) lake, 
(«) terrestrial, and {Hi) wetlands that are seasonally flooded in 
some years, depending on the lake’s water level. 

GBY is located in the southern Hula Valley and assigned to 
the Lower-Middle Pleistocene (marine isotope stages) (18-20), 
~780,000 y ago (12). Excavations and deep core drillings revealed a 
thick sedimentary sequence of the paleo-Lake Hula margin (19) 
deformed by later intensive tectonics (20) and including 26 archae- 
ological layers, among which 15, estimated to represent occupations 
of about 50 ka, are rich (12, 19). Stone artifacts, fossil animal bones, 
and well-preserved organic remains provide data on the Acheulian 
paleoenvironment, ecology, habitat, and cultural realm. 

Reconstructions of the paleolake margin and the flora of its 
diverse adjacent habitats are based on paleobotanical remains 
(12-14, 18) compared with the current local flora (17), while 
taking into account differences imposed by recent anthropogenic 
environmental changes and agricultural activity. Most tree species 
found at GBY (13) still grow today within a 1-km radius from 
GBY, an indication of a Mediterranean climate 780,000 y ago. 

Results 

The botanical remains discussed here have two different sources: 
archaeological layers and geological layers devoid of archaeological 
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remains ( Methods and Table S2). The strata must have been 
rapidly sealed, because in the Mediterranean climate, deterio- 
ration of uncharred plant material exposed to atmospheric 
conditions is swift. 

We consider plants species to be food plants when they are 
consumed by recent rural societies, a minimal criterion because 
the more hardy Paleolithic hominins probably consumed addi- 
tional plant taxa that are not used today. We assume that most of 
the food plant remains were brought to the site deliberately by 
hominins rather than by natural agents. However, an unknown 
proportion of the food plant remains may have arrived without 
hominin intervention. We used food plant frequency in archae- 
ological vs. geological layers as evidence for their deliberate 
collection by the Acheulian inhabitants. Although some plant 
species could have been used for other purposes (fibers, medi- 
cine, fish poisoning, and tool making), we focused on food plants. 
We paid special attention to food plant taxa that appeared in 
conspicuously high proportions in at least four archaeological 
layers and were not common in geological ones. 

Macrobotanical remains were discovered in both archaeological 
and geological layers through the entire 34-m depositional sequence, 
except for a single lignite layer (11-16). Over 100,000 macrobotanical 
fragments were studied. Of these fragments, the minimal number of 


seeds or fruits is 22,714, of which 20,912 were identified to the 
species/genus level. Poor preservation made 1,802 seeds and fruits 
unidentifiable. The identified specimens comprise 117 taxa (78 spe- 
cies and 39 genera), including 48 nonfood taxa (Table S3). 

The food plant assemblage comprises 9,148 plant remains (Fig. 1) 
belonging to at least 55 species. The exact number of species is 
uncertain, because 11 taxa were identified only to the genus level, 
and some of their remains belong to several different but un- 
identified species within these genera. They include nuts, species 
producing USOs, fruits, seeds, and vegetables (Table S4). In 11 
edible species, several organs are eaten (Table S4). Some of the 
archaeological horizons are significantly richer and more diverse in 
edible taxa than others (Fig. 2, Dataset SI, and Tables S5 and S6). 

An obvious question is whether the proportion of food plant taxa 
differs significantly between archaeological and geological layers. Our 
statistical analysis [multidimensional scaling (MDS)] (Methods, Sta- 
tistical Analysis) showed that most horizons are fairly similar with 
respect to plant taxa composition (Fig. SI); only layers 1-4 and V-6 
were outside the boundaries of the 95% ellipse, implying a signifi- 
cantly different plant composition in these two layers relative to all 
others. This MDS model fitted well the observed distance matrix 
(99% of the dispersion was accounted for, and the stress value 
was 0.0098). 



Fig. 1. Food plant remains from GBY. (/\) Quercus sp. f young cupule (layer 11-6 LI); (£) T. natans, upper tip of nut (layer MI-7); (C) Nuphar luteum, seed (layer 11-7); 
(D) B. umbellatus, seed (layer 111-4); (£) Scirpus lacustris, seed (layer 11-9); and (F) Vitis sylvestris, pip (layer IV-7). 
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Fig. 2. Count, frequency, types, and seasonality of food plants found at 
GBY arranged from youngest layer to oldest layer (left to right, respectively) 
(Table S5). (A) Number of taxa and plant foods in the archaeological layers. 
(B) Frequency of food organ types in the archaeological layers. (C) Frequency 
of edible organ types according to seasonality in the richest archaeological 
layers. 


Although the probability of finding food plant remains did not 
differ between archaeological and geological layers when the fre- 
quency of total items or number of plant taxa was considered [like- 
lihood ratio (LR) < 0.7, P > 0.389], the probability of finding items 
relating to key food plants (staples) [Methods, Key (Staple) Food 
Plants \ was significantly greater in archaeological layers than in geo- 
logical layers (LR xj = 5.1, P = 0.024) (Fig. 3). In other words, the 
mean probability of finding key food plants was an order of magni- 
tude higher in archaeological layers than in nonarchaeological 
layers (Fig. 3). 


Discussion 

There is extensive variability in the frequency of food plant taxa and 
organs between the different archaeological layers (Fig. 2, Dataset 
SI, and Tables S5 and S6): four archaeological assemblages (layers 
V-6, II-5, H-6 level 1, and II-6 level 7) are richer than others in plant 
foods (Fig. 2, Dataset SI, and Table S5). This high variability be- 
tween the archaeological layers is not related to the sedimentary 
environment (Methods, Sampling and Sorting of Plant Remains), and 
after calculation of the number of potential types of edible organs, 
vegetables were seen to be the most frequent followed by USOs 
(Fig. 2 and Table S5). However, the actual number of plant remains 
shows that the most nutritious nuts and USO-producing species 
were the most common followed by fruits (Table S7), explaining 
their abundance in layer II-5 and layer H-6 level 1. 

The GBY plant foods include a high diversity of plant organs 
that could have furnished hominins with rich year-round nutri- 
tion. Although the edible plant taxa found are only about 20% of 
the current Upper Jordan Valley food plant taxa (Table SI), they 
include six locally extinct species that existed there during the 
Early-Middle Pleistocene (SI Text, Extinct Species) (15). Of these 
species, Euryale ferox, Sagittaria sagittifolia, and Trapa natans 
could have been used as staple foods, because they are known to 
grow in dense patches in shallow water and are extant crops in 
East Asia (refs. 21, 22, p. 518, and 23). 

The nutritional values of the food plant organs show that nuts 
were the most efficient food source at GBY. Two of them, 
Gorgon nuts (E. ferox) and water chestnut (T. natans), are highly 
nutritious and rich in starch and proteins. Popped seeds of 
Gorgon nuts contain 77% (wt/wt) carbohydrates, 9.7% (wt/wt) 
protein, and 0.1% fat (per 100 g) (24). The seeds of water 
chestnut are composed of 52% (dry wt/dry wt) starch, 15% (wt/wt) 
protein, and 7.5% (wt/wt) fat (23). When oak acorns are added 
as a starch source, combined with Olea fruits and Silybum seeds 
as oil sources, a picture of diverse plant-based nutrition emerges, 
enriched by dozens of other nonstaple food plant taxa. 



Fig. 3. The probability (±95% confidence intervals) of detecting the remains of 
a key food plant in archaeological and geological layers. This probability was 
significantly different between archaeological and geological layers (LR x] = 5.1, 
P = 0.024). 
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Conclusions 

Hominin consumption of plant foods necessitated knowledge 
of their growing localities, seasonality, toxicity, and availability 
(25). Plant phenology in the Hula Valley is mainly dictated by the 
Mediterranean climate (hot, dry summers but rainy winters). The 
seasonal or perennial lake water-level fluctuations either flooded 
or exposed significant areas, causing plants to respond with 
either growth or lack of growth. Winter and spring are rich in 
green vegetables (e.g., Beta vulgaris and Malva nicaeensis), 
whereas in late spring/early summer, many edible seeds and 
fruits are especially available (e.g., Silybum marianum and 
Ziziphus spina-christi). 

Nut and fruit availability is greater in summer and autumn in 
the contemporary Upper Jordan Valley. When the data of all of 
the archaeological layers at GBY are combined, the highest pro- 
ductivity and availability of food plant species are in spring (32 
species) and summer (32 species) (Fig. 2 and Table S5). However, 
each of the archaeological layers contains food plant species of all 
seasons, indicating year-round occupations (SI Text, Seasonality). 

Two factors influenced seasonality and accessibility to food plants 
at GBY. First, there is a considerable altitudinal difference between 
the Hula Valley, the eastern Galilee Mountains, and the Golan 
Heights (a maximal elevation of 800 m over a minimal distance of 
2.5 km), resulting in a difference of several weeks between ripening 
of the same species in the valley vs. the mountains (26). Second, 
moderate lakeshore slopes and oscillations of water level may have 
prolonged or shortened the seasonal availability of aquatic and 
bank food plants, prolonging, for instance, the summer cropping 
season of Sciipus rhizomes, which start growing when water re- 
treats. The duration of palatability differs among species: some have 
a short season, whereas hard nuts have a long one. Acorns, for 
instance, are available for several months, but their seasonality is 
also affected by competition with animals (e.g., wild boar, rodents, 
and birds) that consume large amounts of acorns (27). 

The GBY inhabitants had access to 14 USO-producing species 
(Table S4), many of which could be eaten year round. However, 
the palatability of USOs changes, with highly palatable young 
storage roots/rhizomes gradually becoming fibrous/woody and 
hence, less palatable. Low palatability of USOs can, in certain 
taxa, be handled by roasting/grinding/leaching to deactivate toxins 
or extract starch from fibrous tissues. There is ample evidence for 
the important role of fire at GBY, with its control and repeated 
use shown by burned lithics and charred wood, bark, grains, and 
fruits (28, 29). Fire was instrumental in enriching the value of the 
diet and adding resources that are unpalatable without heating 
(30). Roasting could have enabled the addition to the diet of 
important plant foods: Nuphar lutea rhizomes, Butomus umbella- 
tus, and Sparganium erectum rhizomes (Table SI). Roasting im- 
proved the taste/digestibility of foods, such as oak acorns, 
rhizomes, and young Phragmites australis shoots, common plants 
that could provide considerable amounts of staple foods. Roasting 
also enabled popping of E. ferox (31), and fire could prolong the 
palatability season of tubers of Cyperus rotundas and Scirpus 
maritimus, which are eaten either raw or roasted (32). 

The macrobotanical assemblage of GBY has a surprisingly rich 
vegetal dietary potential that includes dozens of edible species. It 
has provided unprecedented data on Acheulian plant use, illus- 
trating the diversity of the exploited habitats that provided a large 
quantity and wide variety of edible plant organs found in associ- 
ation with a rich material culture. The seasonality of the edible 
plants clearly indicates that the GBY hominins could have occu- 
pied the site year round (Fig. 2 and Table S5). 

We have considered a variety of possible hypotheses to explain 
the abundance of edible plants in the archaeological layers at 
GBY. Among these hypotheses were issues of taphonomy, avail- 
ability, and the possible impact of paleoclimatic changes. Based on 
the results of this study and our past multidisciplinary studies on 


this site (11, 19, 28, 29, 33-38), we suggest that the abundance of 
food plant remains in the archaeological layers is a result of de- 
liberate hominin behavior (Dataset SI). The remains of the key 
food plants are 10 times more abundant in the archaeological 
layers than in the geological ones (Fig. 3). This remarkable dif- 
ference is evident despite the extent of interarchaeological layer 
variability, which also characterizes other components of the ar- 
chaeological horizons (expressed by extensive differences in the 
frequencies of lithic artifacts, fauna, and flora) (11, 29, 39). 

The food plant remains were part of a much more diverse diet 
that probably included additional food plant species, fish (33), 
amphibians, reptiles, birds, mammals, such as fallow deer (40), 
elephants (34), and various aquatic and terrestrial invertebrates. 
The remains of many food plants and animals testily to an ex- 
tremely broad spectrum diet that characterized the site’s many 
archaeological layers. The GBY hominins applied an African 
typotechnology to local raw materials to make lithic artifacts 
that illustrate their African cultural origin but in terms of diet, 
adapted to plant and animal foods that only partly overlap with 
the African ones. The continuous presence of the lake and its 
fluctuating margin habitats together with the surrounding ter- 
restrial habitats provided the GBY hominins with a wealth of 
food resources, enabling repeated occupations of the same locale 
during tens of millennia. These prolonged occupations enabled 
cultural adaptations to the eastern Mediterranean environment. 
The exploitation of new non-African animal and plant taxa tes- 
tifies to Acheulian cognitive/cultural flexibility and adaptability. 
We propose that a wide spectrum of food plants was a perma- 
nent aspect of the preagricultural hominin economy and that its 
infrequent manifestation in archaeobotanical finds probably re- 
flects taphonomic rather than cultural issues. 

Located in the Levantine Corridor and reflecting a Mediter- 
ranean environment and ecology, the food plants of GBY pro- 
vide insight into the resources that enabled hominins’ adaptation 
and survival beyond African habitats, shedding light on the 
mechanisms that enabled their further dispersal into Eurasia. 

Methods 

Plant Context, Taphonomy, and Sampling. Various multidisciplinary studies of 
the record of GBY have shown that the strata were sealed rapidly and hence, 
preserved the components (remains of hominin activity and others) and their 
original spatial organization. These studies include, among others, taphonomic 
studies of malacology (35), crustaceans (41), fish (29), and mammals (36). In 
addition, a wealth of taphonomic information originating from the analyses 
of stone artifacts provides additional insight into the integrity of the assem- 
blages that form the content of the archaeological horizons. It shows that 
there was no winnowing and hence, microartifacts retained their original 
spatial locations; that burned microartifacts were found clustered and thus, 
preserved the original location of hearths; and that some bones from several 
archaeological horizons were successfully conjoined, providing additional 
proof of their taphonomic integrity (29, 36, 42, 43). The presence of perishable 
organic material, such as wood, bark, fruits, and seeds (11, 13, 31), throughout 
the depositional sequence of GBY further testifies to rapid sealing of the 
horizons and hence, the integrity and depositional originality of the edible 
botanical remains. 

Sampling and sorting of plant remains. The excavation area and associated 
geological trenches are located on the left bank of the Jordan River. The 
trenches were dug into the Early and Middle Pleistocene deposits perpendic- 
ularly to the strike of the layers, and hence, their cross-sections provided a 
precise record of their thickness. The Early-Middle Pleistocene deposits were 
not initially visible on the surface of the bank selected for the locations of 
trenching or excavation. 

The trenching revealed a 34-m-thick depositional sequence built of layers 
both containing archaeological remains (archaeological horizons) and devoid 
of them (geological layers) (SI Text, Archaeological and Geological Layers). 
Thus, the discovery of the archaeological horizons was fortuitous and 
moreover, revealed exceptional concentrations of finds (lithics, fossil bones, 
and organic materials). The interlayer differences in frequency of finds are 
indicative of where in the site (and the paleolandscape) the trench or ex- 
cavation was placed — in the center or on the margins. Considering the 
above information, it is possible that layers devoid of stone lithic artifacts 
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and consequently, designated geological may, in fact, represent the edges 
of archaeological layers; the case of a single microartifact (smaller than 2 cm) 
that was found adjacent to a wooden log in layer 11-6 level 14 is an example. 
Theoretically, enlarging the extent of the excavation could have resulted in 
the discovery of finds in the layers presently designated geological. It is 
evident from the stone, bone, and plant finds that the GBY hominins op- 
erated beyond the areas in which stone artifacts and bones were found. 
Sampling of geological layers (devoid of archaeological finds). Samples of 0.5-5 kg 
sediment were obtained from different layers that were exposed in the walls 
of the geological trenches. The samples were placed in sealed nylon bags to 
keep them wet. In the laboratory, the sediments were divided into four 
fractions (0.3-1, 1-2, 2-4, and 4-10 mm) by wet sieving. Items larger than 
1 cm were separated by hand or with large tweezers, and those smaller than 
0.3 mm were lost through the lower sieve. The botanical remains from each 
fraction were separated and sorted by spreading the sediment on trays with 
water and picking them up individually with soft tweezers under a stereo- 
scopic (binocular) microscope at a magnification of up to 25x. 

Sampling of archaeological layers. The entire volume of sediment excavated 
from the archaeological horizons was wet-sieved during fieldwork by a 2-mm 
sieve, and hence, the remains are limited to items larger than 2 mm. The wet- 
sieved sediments were then dried and bagged with their recorded in- 
formation and transported to the Institute of Archaeology for additional 
analysis. Sorting of the sieved sediments yielded rich and varied assemblages, 
such as fruits, seeds, grains, mammalian bones and teeth, fish bones, crab 
skeletons, and specks of charcoal. Many of the seeds and fruits studied here 
(4,199 of 25,835) were retrieved by this procedure. The small-seeded species 
(e.g., Alisma lanceulatum, Chenopodium sp., and Lycopus europaeus) are 
underrepresented in these samples, because they were retrieved only when 
they were stuck or buried in large (>2 mm) lumps of mud. Because the wet- 
sieved sediments were transported from the field with their recorded loca- 
tion, these seeds and fruits could be located within the sediment with a 
precision of 0.5 x 0.5 x 0.5 m. 

Photography. Seeds were photographed to add a visual illustration, serve as a 
basis for future comparison of ancient Levantine flora, and rarely, obtain 
greater confidence in the identifications. Special emphasis was placed on 
seeds of exotic species. Photography was carried out with a scanning electron 
microscope (JEOL model JMS-840) of 10x-100,000x magnification and a 
stereoscopic microscope (Olympus model SZX12) of up to 90x magnification. 

Seeds were cleaned by immersion in water using paintbrushes and needles 
to prepare them for photography. Seeds prepared for SEM photography 
were pasted on a stab and coated with gold for 10-20 min (depending on 
their size, shape, and texture). Waterlogged seeds are difficult to dry 
without destruction of shape and therefore, were fixed by substituting the 
water with organic materials (critical point drying method). SEM photog- 
raphy was performed at the Faculty of Life Sciences, Bar-llan University with 
the help of Yakov Langsam. The microscope digital photographs were 
processed by image-editing software (Paint Shop Pro-7). In cases of speci- 
mens larger than 1 cm (the maximal size of the SEM chamber used), two 
parts of the seed were photographed successively and later combined. In 
cases where it was impossible to achieve the same focus for the two pictures, 
the determination of the size of the object was slightly affected. In such 
cases, a small space was left between the two photographs. Tiny pieces of 
carbon glue strips were attached to the side of the specimen or filled spaces 
between specimens to overcome difficulties in SEM photography resulting 
from the height of large specimens. 

Taxonomic Identification. 

Seed and fruit identification and calculation of the number of specimens. Water- 
logged seeds are extremely sensitive to dryness and can lose their shape easily. 
Consequently, their processing requires gentle handling during identification 
and preparation for photography. Furthermore, many of the seeds were found 
broken, making them even more difficult to identify. The identification process 
relies on experience and familiarity with the morphology and anatomy of the 
seeds, fruits, and other plant parts of the local flora. This experience is based on 
the study of the morphological characteristics of the seeds known for the 
taxonomic group or groups. Identification within the groups is based on ex- 
amination of the seeds and fruits with the aid of a reference collection of the 
plants of Israel and the Middle East (Faculty of Life Sciences, Bar-llan University) 
and publications including illustrated atlases of the plants of Israel and flora of 
the region (44-48) and worldwide (49-51) and was complemented at times by 
target-oriented field work in the vicinity of the site. 

1. Ungar PS, Grine FE, Teaford MF (2006) Diet in early Homo: A review of the evidence 

and a new model of adaptive versatility. Annu Rev Anthropol 35:209-228. 


The taxa classified as vegetable-producing species were identified by their 
seeds. Therefore and because in some of them, more than one organ can be 
used as food, we did not have a direct method for counting the number of 
specimens for statistics. 

Plant remains that were identified only to the family or tribe level were 
defined as unidentified and were not considered to be a component of the 
number of taxa, which is limited to remains identified to the genus or species 
level. When two remains were identified, one, because of its bad preser- 
vation, to the genus level and the other to a specific species of the same 
genus, they were considered as two different taxa. When two known variants 
of a species were identified, we considered them as one taxon. Aegilops 
geniculata/peregrina and Aegilops cf geniculate as well as Ziziphus lotus/ 
spina-christi and Z. spina-christi are pairs of related species that provide the 
same edible organ in the same season. Therefore, when counting the 
number of edible plants, we considered each pair as one taxon. 

Here, we consider leaves and young shoots as vegetables, nuts as either 
very large seeds with a hard coat or fruits with a hard shell, and Gramineae 
seeds as grains. We determine seasonality according to the phenology of 
plants under Mediterranean conditions (some plant foods can be gathered in 
more than one season). 

Key (staple) food plants. Key food plants include the following taxa: E. ferox, 
Quercus acorns, 5. marianum, Olea europea, T. natans, S. sagittifolia, and 
Typha. These taxa were selected because of their high nutritional value and 
the possibility of easily gathering large amounts of their edible organs. 

Statistical Analysis. We used MDS (PROXCAL) (52) to cluster layers based on 
the relative abundance of plant remains. MDS is a robust approach for vi- 
sualizing the pattern of proximities (i.e., similarities or distances) among a 
set of objects. To accommodate for differences in sediment volume exam- 
ined from each layer, we calculated the proportion of items recovered from 
each plant species from the total items recovered in each layer. The relative 
occurrence of the remains of 36 identified plant species was used for cal- 
culating the Chebychev distance between each pair of layers. Stress (i.e., the 
degree of correspondence between the distances among points calculated 
by the MDS map and the input matrix) and Shepard diagram (i.e., a scat- 
terplot of the input against the output proximities) were used as the mea- 
sure of fit between the observed and calculated distance matrices. To cluster 
layers into groups, we used 95% confidence ellipse. 

We used plant edibility (i.e., edible or not for humans) as a binary de- 
pendent variable in a logistic regression under the framework of generalized 
linear models. Model distribution was set as binomial, and the link function 
was set as logit. Each pair of rows in our data stored the identified plant 
information for one layer: the top for food plants and the next for nonfood 
plants. Soil type, classification of the layer (i.e., archaeological or geological), 
and chronological order of the layer were used as independent variables in 
our models. The frequency of food items in each layer was corrected for the 
volume of sediment examined from each layer. LR y 1 2 was used for evaluation 
of the effect of the above variables on the probability of finding plant food 
items. Calculations were performed using SPSS (version 22; SPSS Inc.) and 
JMP (version 12; SAS Inc.). 

The type of sediment had a significant effect on the probability of re- 
covering food plant remains (all plant items: LR y\ = 398.1, P < 0.001; item of 
key food plants: LR y\ = 49.2, P < 0.001). The probability of recovering food 
plant remains was significantly higher (0.54-0.64) in the storm beach and BC 
(a contact between black mud and coquina) sediments compared with all 
other sediment types (0.01-0.03). However, this trend was observed only in 
the geological layers (all plant items: LR y% = 396.4, P < 0.001; items of key 
food plants: LR y% = 44.1, P < 0.001). In archaeological layers, the probability 
of recovering food plant remains was independent of soil type (all plant 
items: LR y 2 = 1 .2, P = 0.757; items of key food plants: LR yl = 1 .4, P = 0.737). 
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